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An  efficient  solution  to  reducing  NOx  formation  is  to  maintain  a  relatively  low  flame  temperature.  This 
can  be  achieved  by  mixing  reactants,  prior  to  combustion,  with  chemically  inert  diluents  such  as  cooled 
combustion  products.  In  such  diluted  combustion  systems,  the  flame  temperature  decreases  because  of 
thermal  ballast,  limiting  NOx  production.  This  work  focuses  on  modeling  the  specifics  of  this  combustion 
regime  in  confined  combustors.  To  characterize  the  dilution  of  reactants  by  burnt  gases,  the  importance 
of  complex  chemistry  effects  is  emphasized  and  taken  into  account  using  a  detailed  chemistry  tabulation 
approach.  This  approach  extends  the  flamelet/progress  variable  formulation  by  including  information 
about  the  intensity  of  internal  dilution  rates  and  heat  losses.  A  turbulent  combustion  model  is  then  devel¬ 
oped  in  a  large  eddy  simulation  (LES)  framework.  The  combustion  model  is  validated  by  considering  two 
combustor  configurations,  namely  an  adiabatic  burner  and  a  combustor  having  isothermal  walls  -  both 
operating  under  highly  diluted  combustion  conditions.  Simulation  results  are  in  good  agreement  with 
experimental  data,  confirming  the  importance  of  detailed  chemistry  information  and  the  validity  of 
the  tabulation  approach  to  LES  application  to  diluted  combustion. 

©  2014  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 


1.  Introduction 

The  energetic  efficiency  of  combustion  systems  can  be 
enhanced  by  transferring  heat  from  exhaust  products  to  the  fresh 
gases  by  means  of  regenerative  heating.  Unfortunately,  the  result¬ 
ing  increase  in  the  reactant  temperature  has  an  adverse  effect  on 
the  formation  of  nitrogen  oxides  (NOx).  An  efficient  solution  to 
reducing  NOx  formation  is  to  maintain  a  relatively  low  flame 
temperature.  This  can  be  achieved  by  mixing  reactants  with  chem¬ 
ically  inert  diluents  such  as  combustion  products.  In  “diluted  com¬ 
bustion”  systems,  the  flame  temperature  decreases  because  of 
thermal  ballast,  limiting  NOx  production.  These  systems  promote 
different  combustion  regimes  such  as  MILD  (moderate  or  intense 
low-oxygen  dilution)  combustion  [1,2],  fiameless  oxidation  [3] 
and  high-temperature  air  combustion  (HiTAC)  [4,5], 

In  diluted  combustion  technologies,  reactants  are  diluted  with 
large  amounts  of  burnt  reaction  products  prior  to  combustion, 
which  enables  flame  stabilization  under  lean  conditions,  thereby 
avoiding  high-temperature  regions  that  promote  enhanced 
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thermal  NOx  formation.  Experimental  and  numerical  studies  have 
been  conducted  in  confined  systems  such  as  HiTAC  combustion 
chambers  [6,7],  the  IFRF  (International  Flame  Research  Founda¬ 
tion)  semi-industrial-scale  configurations  [8-10],  and  reverse  flow 
configurations  [11,12].  These  diluted  combustion  technologies 
exhibits  strong  recirculation  zones  that  enhance  the  mixing 
between  fresh  and  burnt  streams.  Local  stoichiometric  conditions 
are  avoided,  so  that  NOx  production  is  dramatically  decreased.  It 
is  noteworthy  that  diluted  combustion  can  be  achieved  without 
internal  recirculation  of  burnt  gases:  for  example,  Dally  et  al. 
[13[  reproduced  oxygen-diluted  and  MILD  combustion  regimes  in 
a  Jet  in  Hot  Coflow  (JHC)  experiment. 

A  schematic  of  the  model  problem  of  diluted  combustion  with 
internal  product  gas  recirculation  in  a  confined  burner  geometry 
is  shown  in  Fig.  1.  The  schematic  illustrates  the  modification  in 
the  composition  of  reactants,  which  strongly  varies  with  recircula¬ 
tion  and  dilution  by  burnt  gases. 

A  major  issue  in  the  modeling  of  diluted  combustion  is  the 
pronounced  sensitivity  of  the  flame  structure  to  the  reaction  chem¬ 
istry  [14,15].  To  illustrate  the  importance  of  the  detailed  chemistry 
when  a  gas  mixture  is  subjected  to  dilution  by  hot  reaction 
products,  a  series  of  adiabatic  laminar  flames  are  computed.  These 
simulations  consider  chemical  representations  of  increasing 
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Nomenclature 

</> 

equivalence  ratio 

Latin 

<t>c 

global  equivalence  ratio 

C 

normalized  progress  variable 

X 

scalar  dissipation 

D 

molecular  diffusivity 

* 

vector  of  chemical  quantities 

H 

specific  mixture  enthalpy 

6) 

vector  of  chemical  reaction  rates,  expressed  in  s_1 

Sz 

mixture  fraction  unmixedness  factor 

cok 

chemical  reaction  rate  of  species  k,  expressed  in  s_1 

T 

temperature 

Yc 

progress  variable 

Superscripts 

Yd 

dilution  variable 

Dil 

diluent  stream 

Ye 

elemental  mass  fraction 

F 

fuel  stream 

Yk 

species  mass  fraction  of  species  k 

Ox 

oxidizer  stream 

Z 

mixture  fraction 

0 

undiluted  conditions 

Z 

mixture  fraction  normalized  with  respect  to  undiluted 

compositions 

Subscripts 

G 

global  equivalence  ratio  conditioning 

Greek 

la=0 

undiluted-conditioned  quantity 

(X 

dilution  parameter 

I/M) 

quantity  evaluated  without  heat  losses  (/J  =  0) 

p 

heat  loss  parameter 

1/1=1 

quantity  evaluated  for  maximal  heat  losses  (/?=!) 

A 

progress  parameter 

P 

density 

complexity,  namely  the  infinitely  fast  chemistry  approximation,  a 
one-step  irreversible  chemistry,  and  a  detailed  reaction  sequence. 
The  problem  configuration  consists  of  a  steady-state  counterflow 
methane/air  diffusion  flame,  and  different  values  of  strain  rates 
(including  the  stable  and  the  unstable  branch)  are  considered.  Both 
reactant  streams  are  diluted  with  burnt  products,  having  a  compo¬ 
sition  that  is  identical  to  that  of  the  MILD  combustor  configuration 
studied  by  Verissimo  et  al.  [16]  (see  Section  5.1).  Chemical  trajec¬ 
tories  in  C02— CH4  composition  space  are  presented  in  Fig.  2  for  dif¬ 
ferent  dilution  levels.  Results  obtained  with  the  detailed  GRI  3.0 
mechanism  [17]  are  shown  by  solid  lines,  the  dashed  line  corre¬ 
sponds  to  the  one-step  irreversible  chemistry,  and  the  limiting  case 
of  infinitely  fast  chemistry  is  shown  by  the  symbol.  From  this  fig¬ 
ure  it  can  be  seen  that  the  detailed  chemistry  solution  exhibits  sev¬ 
eral  possible  trajectories,  which  depend  on  the  level  of  dilution,  in 
contrast,  the  single-step  chemistry  trajectory  is  not  sensitive  to  the 
dilution  level,  so  the  CH4  mass  fraction  exhibits  a  linear  depen¬ 
dence  on  Y Co2.  Since  the  infinitely  fast  chemistry  model  assumes 
that  the  mixture  is  always  at  equilibrium  state,  the  chemical  trajec¬ 
tory  reduces  to  a  single  point  for  a  given  value  of  equivalence  ratio. 
Results  from  the  detailed  chemistry  solution  show  that  the  chem¬ 
istry  is  affected  by  dilution,  impacting  fundamental  flame  proper¬ 
ties,  including  flame  structure,  species  composition,  and  pollutant 
emission. 

An  attractive  strategy  for  including  detailed  chemistry  effects 
using  moderate  CPU  resources  are  tabulated  chemistry  techniques 
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Fig.  1.  Model 

problem 

of  a  diluted  combustion 

configuration  with  internal 

recirculation  of  burnt  gases.  Considering  the  idealized  problem  of  methane/oxygen 
combustion,  the  initially  separated  reactants  mix  and  react  to  produce  C02  and 
H20.  Reaction  products  are  then  recirculated  where  they  modify  the  composition  in 
the  oxidizer  and  fuel  streams.  The  product  gas  dilution  level  is  denoted  by  a. 


[18-22].  Among  these,  the  flamelet  model  for  nonpremixed  com¬ 
bustion  assumes  that  a  turbulent  flame  can  be  decomposed  into 
a  collection  of  one-dimensional  flame  elements  [18],  Each  flamelet 
is  then  represented  by  a  reaction-diffusion  element  that  is  con¬ 
structed  between  oxidizer  and  fuel  streams.  In  the  original  formu¬ 
lation,  proposed  by  Peters  [18],  the  fuel  and  oxidizer  composition 
are  assumed  to  be  constant  for  each  flame  element.  By  construc¬ 
tion,  this  two-stream  formulation  is  not  able  to  account  for  effects 
of  reactant  dilution  by  burnt  gases  on  the  chemical  flame  structure. 
To  overcome  this  issue,  a  three-stream  flamelet-progress  variable 
(FPV)  approach  has  recently  been  developed  [23],  This  model 
was  applied  to  large  eddy  simulations  of  a  Jet-in-Hot-Coflow 
(JHC)  burner  [23,24],  in  which  the  burner  was  operated  in  the 
recirculation-free  adiabatic  MILD  operating  regime.  The  dilution 


Fig.  2.  Chemical  trajectories  in  C02-CH4  state  space  for  different  dilution  levels, 
illustrating  the  sensitivity  of  the  fuel  conversion  to  the  reaction  chemistry  and  the 
dilution.  Here,  the  diluent  is  composed  of  the  equilibrium  product-gas  composition 
for  an  equivalence  ratio  of  0  =  0.58.  Trajectories  are  extracted  from  steady  laminar 
counterflow  diffusion  flame  computations,  from  the  pure  mixing  line  to  the  fully 
burnt  states. 
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is  provided  by  a  hot  coflow  stream,  so  the  diluted  reactant  compo¬ 
sition  does  not  vary  in  the  combustion  chamber.  A  model  that  uses 
two  independent  mixture  fractions  is  therefore  sufficient  to  cap¬ 
ture  the  flame  properties  [23],  Nonetheless,  when  the  dilution  pro¬ 
cess  is  induced  by  strong  internal  recirculation,  the  local 
composition  of  the  diluted  fuel  or  oxidizer  stream  changes  in  time 
and  space,  and  the  use  of  passive  scalars  to  describe  the  mixture 
variations  is  not  adequate  any  more. 

The  present  work  presents  a  modeling  route  to  account  for  the 
effect  of  fresh-gas  dilution  by  burnt  products  on  the  dynamics  and 
structure  of  turbulent  flames.  This  approach  is  applied  to  the  pre¬ 
diction  of  “diluted  combustion”  systems,  including  for  instance 
MILD  (moderate  or  intense  low-oxygen  dilution)  combustion 
[1,2  or  flameless  oxidation  [3]  regimes.  To  this  end,  a  tabulation 
methodology  based  on  the  FPV  method,  which  accounts  for  effects 
of  internal  dilution  on  the  chemical  flame  structure,  is  first  intro¬ 
duced.  Since  heat  losses  are  important  in  most  practical  combus¬ 
tors,  the  impact  of  these  effects  on  the  chemistry  will  also  be 
considered.  Two  tabulation  strategies  that  differ  in  the  representa¬ 
tion  of  heat  transfer  and  dilution  are  considered.  A  sub-grid-scale 
model  for  the  turbulence/flame  interaction  is  proposed,  and  a 
semiempirical  model  is  developed  to  represent  the  production 
and  internal  recirculation  of  diluents.  The  capability  of  this  LES 
combustion  model  is  demonstrated  in  two  combustor  configura¬ 
tions,  namely  an  adiabatic  burner  and  a  combustor  that  utilized 
isothermal  wall  treatment.  These  configurations  were  experimen¬ 
tally  investigated  at  the  Institute  Superior  Tecnico  of  Lisbon 
[25,16],  Both  burners  feature  severe  internal  recirculation,  very 
low  pollutant  emissions,  and  high  combustion  efficiency.  Compar¬ 
isons  with  experimental  data  are  then  carried  out,  and  the  validity 
of  the  proposed  model  is  discussed. 

2.  Modeling  of  diluted  combustion 


yco2 

Fig.  3.  Yco  evolution  as  a  function  of  YCo2  during  autoignition  of  a  PSR  for  several 
residence  times  z.  The  solid  trajectory  corresponds  to  the  0-D  reactor  computation 
(with  t  — *  oo  in  Eq.  (1)). 

corresponds  to  a  given  residence  time  r  of  the  gases  inside  the  reac¬ 
tor.  The  solid  line  represents  the  chemical  trajectory  of  a  homoge¬ 
neous  reactor  (solution  of  Eq.  (1)  with  r  — >  oo).  For  a  long 
residence  time  (t  =  100  ms),  the  chemical  trajectory  rapidly  con¬ 
verges  toward  the  closed  reactor  solution,  and  the  fresh  gases  dilu¬ 
tion  by  burnt  gases  does  not  significantly  affect  the  CO  production. 
However,  for  shorter  residence  times  (t  =  10  ms  and  t  =  1  ms),  tra¬ 
jectories  differs.  The  dimension  of  the  attractive  manifold  increases 
and  a  chemistry  tabulation  function  limited  to  a  single  progress  var¬ 
iable  (here  C02)  is  not  sufficient  to  capture  the  chemical  processes. 
A  solution  to  include  dilution  effects  in  the  FPV  model  is  presented 
in  the  following  section. 


2.1.  Influence  of  fresh  gas  dilution  on  the  chemical  trajectories 


2.2.  Chemistry  tabulation 


The  equilibrium  thermochemical  state  is  identified  from  the 
element  composition  and  the  enthalpy  h.  Under  unity  Lewis  num¬ 
ber  assumption,  the  pair  (Z,H),  where  Z  is  the  mixture  fraction  and 
H  the  mixture  enthalpy,  is  then  sufficient  to  capture  the  equilib¬ 
rium  state  of  the  mixture,  whether  or  not  fresh  gases  have  been  di¬ 
luted  by  burnt  gases.  This  is,  however,  not  the  case  for  reaction 
processes  that  are  strongly  affected  by  the  rate  of  fresh  gas  dilu¬ 
tion.  To  illustrate  this,  we  consider  a  0-D  constant-pressure  perfec- 
tively  stirred  reactor  (PSR),  which  is  described  by  the  following  set 
of  governing  equations: 


f>d0t=r(Y"-Yk)+P°h 

u  \fc-l  k=l  /  k=  1 


0) 


The  inlet  mass  flow  rate  m  and  the  residence  time  t  are  related  to 
the  volume  V  of  the  reactor  through  the  relation  m  =  pV / t.  The 
reactor  is  initially  fed  with  a  mixture  of  fresh  gases  at  temperature 
T(t  =  0)  and  composition  Yk(t  =  0).  Under  constant-pressure  condi¬ 
tions,  this  mixture  reaches  the  equilibrium  state,  denoted  by  super¬ 
script  b.  At  t  =  0,  gases  of  temperature  Tb  and  composition  Ybk  are 
injected  into  the  reactor  and  dilute  the  fresh  gas  mixture.  As  the 
elemental  composition  of  the  mixture  is  kept  constant  during  the 
dilution  process,  the  thermochemical  equilibrium  state  is  not  af¬ 
fected.  A  series  of  PSR  simulations  are  performed  with  the  following 
initial  gas  composition  and  temperature:  Y  o2  =  0.14253, 
Yc h4  =  0.00214.  VN2  =  0.75633,  yH2o  =  0.099,  and  T  =  1330  K,  and 
using  the  GRI  3.0  mechanism  [17],  Figure  3  shows  the  solutions  of 
Eq.  (1)  projected  in  the  (YCo,Yco2)  chemical  subspace.  Each  curve 


2.2.3.  Flamelet/progress  variable  approach 

In  the  FPV  approach,  the  turbulent  flame  is  locally  seen  as  a  col¬ 
lection  of  one-dimensional  laminar  counterflow  diffusion  flames. 
Assuming  steady  combustion  and  neglecting  preferential  diffusion 
effects,  the  structure  of  these  flame  elements  is  then  described  by 
the  solution  of  the  flamelet  equations,  which  are  here  written  as 


dZ2 


(2) 


where  Z  is  the  mixture  fraction,  i/<  is  a  vector  composed  of  all  species 
mass  fractions  Y  and  temperatures  T,  and  the  vector  &>  denotes  their 
respective  source  terms.  =  2D  |  VZ|2  is  the  scalar  dissipation  rate 
of  the  mixture  fraction,  which  introduces  the  molecular  diffusivity 
D  of  the  mixture.  Since  the  focus  of  this  investigation  is  the  diluted 
combustion  of  a  methane/air  mixtures,  we  assume  equal  diffusivi- 
ties  for  all  species. 

To  accommodate  variations  in  the  reactant  stream  composi¬ 
tions,  Eqs.  (2)  are  supplemented  by  the  following  boundary 
conditions: 


+(Z  =  0)  =  <l>°\  (3) 

!>{Z=  1)  =  ^F,  (4) 


where  the  superscripts  F  and  Ox  denote,  respectively,  boundary 
conditions  in  the  fuel  stream  and  in  the  oxidizer  stream.  Thermo¬ 
chemical  quantities  are  then  functions  of  the  mixture  fraction  and 
the  scalar  dissipation  rate.  However,  this  parameterization  does 
not  ensure  a  unique  thermochemical  state  space  representation, 
since  for  given  values  of  xz,  stable,  unstable,  and  unburned  states 
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may  coexist.  To  overcome  this  issue,  a  reaction  progress  parameter 
A  has  been  introduced  [22,26],  which  provides  a  unique  representa¬ 
tion  of  the  entire  thermochemical  state  space.  The  progress  param¬ 
eter  is  defined  through  a  progress  variable  Yc,  which  can  be 
evaluated  as  the  solution  of  a  minimization  problem  [27].  A  defini¬ 
tion  for  Yc  in  the  context  of  methane/air  combustion  is 


flamelet  archetype  that  is  used  for  the  chemistry  tabulation.  The 
counterflow  flamelet  boundary  conditions  now  read 

+(Z=  0)  =  ^0x(a,y5),  (10) 

f(Z=l)  =  ^F(«,/0,  (11) 


yc  =  yco2  +  v,co  +  yH2  +  ni2o.  (5) 

Yc  and  Z  are  by  definition  statistically  dependent.  To  facilitate 
the  estimation  of  the  joint  presumed  probability  density  function 
(PDF)  for  the  turbulent  closure,  the  variable  A  is  introduced.  A  is 
evaluated  by  conditioning  Yc  at  a  given  value  of  the  mixture  frac¬ 
tion  Z", 

A  =  YC|Z.,  (6) 

and  the  thermochemical  vector  ^  is  a  function  of  mixture  fraction 
and  progress  parameter: 


where  a  is  the  dilution  parameter  and  p  is  the  heat  loss  parameter. 

A  schematic  of  representative  flame-configurations  and  defini¬ 
tions  of  the  dilution  and  heat-loss  parameters  is  shown  in  Fig.  4. 

The  dilution  parameter  a  is  defined  as  the  mass  ratio  between 
fresh  reactants  and  diluents:  a.  =  0  in  the  absence  of  diluents  and 
a  =  1  when  the  mixture  is  saturated  with  diluents.1  To  identify  the 
composition  of  the  diluent  stream,  we  introduce  Y^'1,  denoting 
the  mass  fraction  of  species  k  in  the  diluent  stream.  In  practice, 
the  diluent  mixture  composition  Y°‘'  is  evaluated  from  the  thermo¬ 
chemical  equilibrium  state  of  the  mixture  under  global  equivalence 
ratio  conditions.  Denoting  by  Yfk  and  Y°x  the  mass  fractions  of  species 
k  in  the  fuel  and  oxidizer  streams,  the  definition  of  a  reads 


*  =  F+(Z,  A).  (7) 

Since  the  balance  equation  for  A  is  difficult  to  close  for  turbu¬ 
lent  flows,  it  is  preferable  to  solve  a  transport  equation  for  the  pro¬ 
gress  variable  Yc.  The  function  Ty,  ,  which  relates  Yc  to  Z  and  A,  is 
identified  by  using  Eqs.  (5)  and  (7): 

YC  =  ^,(Z,A).  (8) 

Every  thermochemical  quantity  contained  in  ip  evaluated  in  the 
steady  flamelet  equations  is  therefore  expressed  as 


yF  yF,0  yOx  yOx,0 

a  =  _^ _ h _ =  _ Y_k _  (12) 

YDk\p)-Y™  Y^-Y?*'0’ 

and  the  superscript  “,0”  denotes  the  undiluted  mixture. 

The  heat  loss  parameter  p  is  introduced  to  consider  nonadia- 
batic  effects.  This  parameter  is  equal  to  zero  under  adiabatic  con¬ 
ditions  and  [S  —  1  when  heat  losses  are  maximum.  By  denoting  as 
HDl1  the  enthalpy  (sum  of  chemical  and  sensible  enthalpies)  of 
the  diluent,  the  definition  of  p  is  written  as 


^  =  ^(Z^y](Z.Yc))  =  g^Z.Yc). 


(9) 


2.2.2.  Diluted  flamelet/progress  variable  approach 

Capturing  dilution  effects  in  a  flamelet-based  tabulated  chemis¬ 
try  framework  implies  first  modeling  the  mixing  of  the  fresh  reac¬ 
tants  with  the  burnt  product  gases.  As  mentioned  previously, 
assuming  equilibrium  state,  the  burnt  gas  composition  and  tem¬ 
perature  depend  only  on  the  mixture  fraction  and  enthalpy  under 
a  unity  Lewis  number  assumption.  These  reaction  products  may 
mix  at  different  rates  with  the  fresh  oxidizer  and/or  the  fuel 
streams,  multiplying  the  number  of  possible  configurations.  The 
degrees  of  freedom  underlying  these  processes  are  so  numerous 
that  a  direct  tabulation  of  flamelets  for  all  dilution  configurations 
is  in  practice  not  feasible.  However,  reasonable  simplifications 
can  be  made  as  in  confined  burners,  the  dilution  is  accomplished 
by  large-scale  recirculation  of  burnt  gases.  Therefore,  diluted  com¬ 
bustion  can  be  represented  as  a  three-stream  problem,  consisting 
of  a  fuel  stream,  an  oxidizer  stream,  and  an  internal  reaction  prod¬ 
uct  gas  stream.  Because  of  large  residence  time  in  the  recirculation 
zones,  the  burnt  gas  composition  corresponds  to  the  chemical 
equilibrium  of  a  fuel/oxidizer  mixture  that  is  characterized  by 
the  global  equivalence  ratio  </>G.  These  large  time  scales  also  favor 
heat  exchanges  between  recirculating  burnt  gases  and  burner 
side-walls.  Therefore,  the  diluent  composition  is  in  practice  mainly 
a  function  of  the  global  equivalence  ratio  and  the  level  of  wall  heat 
transfer. 

The  study  by  Abtahizadeh  et  al.  [28]  showed  that  the  dilution 
process  has  a  direct  effect  on  the  flame  structure.  In  terms  of  chem¬ 
istry  tabulation,  accounting  for  the  influence  of  diluent  gradients 
considerably  increases  the  dimensionality  of  the  tabulated  mani¬ 
fold.  In  the  approach  developed  here,  we  will  neglect  the  influence 
of  these  gradients  on  the  flame  structure. 

By  these  assumptions,  boundary  conditions  of  the  system  of  Eq. 
(2)  are  reformulated  to  account  for  the  effect  of  dilution  on  the 


rrDil  rrDil  i 

o  h  h  \P=0 

f_rDil  |  rrDil  i 

“  l/?=l  ^  I/M) 

(13) 

From  this  definition,  the  composition  and  enthalpy  of  the 
diluted  reactant  streams  (fuel  and  oxidizer)  can  be  deduced: 

Yj(«,/0  =  a(Y“(/0-Yjfl)+Yjfl, 

(14) 

Hf(oc,  p)  =  a(HDi\p)  -  Hfo^J  +  Hf  0, 

(15) 

Y°x(a,/S)  =  a(Y“'(/?)  -  Y°x0)  +  Y°x0, 

(16) 

H°*(ct,p)  =  x(hd"{P)  -H0x°)  +H0x0, 

(17) 

and  Hm  is  written  as 

H°"(P)  =  p(Hm\M  -  Hd'\= o)  +  HDil|;j=0.  (18) 

The  minimal  enthalpy  of  the  diluent  enthalpy  HDll|^=1  corre¬ 
sponds  to  the  enthalpy  of  the  diluent  mixture  that  is  cooled  to  a 
specified  minimum  temperature  T™n.  In  the  following  this  temper¬ 
ature  is  associated  with  the  burner  wall  temperature. 

When  all  thermochemical  quantities  are  parameterized  in 
terms  of  mixture  fraction,  progress  variable,  dilution  parameter, 
and  heat-loss  parameter,  t p  evolves  in  a  four-dimensional  subspace. 
In  analogy  to  Eq.  (9),  this  is  written  as 


<l>  =  9,/,(Z,  Yc,  a,  p),  (19) 

and  this  formulation  is  referred  to  as  diluted  flamelet/progress  var¬ 
iable  (DFPV)  approach. 


1  Unlike  [23],  this  dilution  parameter  is  not  a  second  mixture  fraction.  It  tracks  the 
deviations,  induced  by  dilution,  of  the  chemical  trajectories  projected  in  a  constant  Z 
plane. 
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Fig.  4.  Schematic  representation  of  the  different  constraints  imposed  on  temperature  and  species  of  the  boundary  conditions  used  to  generate  counterflow  diffusion  flames; 
left:  FPV  model;  right:  DFPV  model. 


2.3.  Relations  between  controlling  scalars 

All  control  parameters  that  appear  in  Eq.  (19)  are  evaluated 
from  transported  quantities.  The  mixture  fraction  Z  is  defined  from 
a  set  of  element  mass  fractions  Ye  [29]: 


in  mixture-fraction  space  decreases  with  increasing  dilution,  until 
it  reaches  a  singular  point  that  is  defined  by  the  pure  diluent 
mixture. 

The  heat  loss  parameter,  p,  is  linked  to  the  enthalpy,  H,  which  is 
expressed  in  the  Z-mixture  fraction  space  as 


-  C(«) 

Ve(«)  -  ye»  ' 


(20) 


Since  the  fuel  and  oxidizer  elemental  species  mass  fractions 
(Yj(a)  and  Y°x(cc)  respectively)  are  functions  of  the  dilution  ratio, 
the  balance  equation  for  Z  exhibits  unclosed  dependencies.  To 
avoid  the  dependency  on  the  dilution  parameter,  we  introduce  a 
mixture  fraction  formulation,  which  is  defined  with  respect  to 
the  undiluted  mixture  composition,  and  therefore  independent  of 
a.  This  quantity  is  denoted  by  Z  and  has  the  following  definition: 


Z  =  - 


Ye  -  Y  „ 


Y  I  _  Y 

1  e  la-0  1e 


(21) 


As  such,  Z  can  be  directly  evaluated  by  solving  a  conserved  sca¬ 
lar  transport  equation.  Furthermore,  by  introducing  the  mixture 
fraction  Zc,  which  is  evaluated  under  global  equivalence  ratio  con¬ 
ditions  (20)  and  (21)  are  related  as 


Z  —  ccZq 
1  -  a 


(22) 


Figure  5  illustrates  temperature  profiles  as  a  function  of  mixture 
fraction  Z  and  different  values  of  reaction  progress  Yc  and  dilution 
levels.  In  this  configuration,  the  oxidizer  consists  of  preheated  air 
at  600  I<  and  the  fuel  stream  consists  of  methane  at  ambient  tem¬ 
perature.  The  diluent  is  composed  of  burnt  reaction  products  at  an 
equivalence  ratio  of  =  0.77  without  any  heat  losses  (p  =  0).  The 
chemistry  is  described  by  the  GRI  3.0  mechanism  [17].  The  flame 
structure  is  obtained  from  steady  laminar  flamelet  computations 
using  the  FlameMaster  code  [30].  This  figure  shows  that  the  range 


H(Z,  a,p)  =  Z  jffF(a,  /?)  -  H0x(a,  /?)J  +H0x(a,j8).  (23) 

Combining  Eqs.  (15),  (17),  and  (22),  and  using  the  definition  of  p 
from  Eq.  (18),  leads  to 


H  -  (Ci o  ~  -  aZc)  -  (1  -  a)H°l0  -  «HDil|,=0 

tx(HD%,-Hm\^0) 


(24) 


Note  that  this  expression  is  valid  for  a.  #  0,  which  is  consistent 
with  the  underlying  assumption  that  heat-loss  and  dilution  effects 
are  coupled.  Consequently,  a  =  0  implies  p  =  0. 

For  this  approach,  the  rate  of  dilution  needs  to  be  estimated 
from  the  evaluation  of  the  reactant  composition  prior  to  combus¬ 
tion.  Dilution  with  burnt  gases  as  well  as  chemical  reactions  will 
affect  the  progress  variable  Yc,  defined  from  Eq.  (5).  To  discrimi¬ 
nate  the  progress  of  reaction  from  dilution  effects,  the  normalized 
form  of  the  progress  variable,  C,  is  introduced, 


c  Yc-Y°c(Z,x,p) 

Yec\Z,uJ)-Y°c(Z,a,py 

where  Y®q  and  Y°  represent,  respectively,  the  value  of  Yc  at  equilib¬ 
rium  and  in  the  unburned  mixture  diluted  at  a  rate  a. 

The  dilution  variable  Yd  is  defined  from  a  combination  of  spe¬ 
cies  mass  fractions,  chosen  such  that  Yjj°  =  0  and  YjX  °  =  0.  The  fol¬ 
lowing  relation  between  Yd  and  a  is  then  obtained  from  Eq.  (12): 


Yd 

YDda(py 


(26) 


Yd‘'(/3)  denotes  the  value  taken  by  the  dilution  variable  Yd  in  pure 
diluent,  i.e.  at  equilibrium  and  for  global  equivalence  ratio  mixing 


Fig.  5.  Temperature  profiles  as  function  of  mixture  fraction  Z  for  different  levels  of  reaction  progress  and  increasing  levels  of  dilution  levels  a  (from  left  to  right). 
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conditions.  In  the  following  simulations  the  dilution  variable  is  de¬ 
fined  as  Yd  =  YCo  +  Yco2-  From  gas  equilibrium  computations,  we 
observed  that  this  combination  ensures  a  very  weak  dependency 
of  Y°'1  on  the  heat-loss  level  p.  The  definition  of  a.  is  then  simplified 
to 


a  = 


yd 

y“'(/?  =  0)' 


(27) 


This  simplification  was  tested  for  the  experimental  configura¬ 
tions  and  operating  conditions  investigated  further,  and  was  found 
adequate  for  all  mixture  compositions.  Yet  this  assumption  should 
be  investigated  further  for  application  to  heavy  hydrocarbon  fuels, 
which  has  not  been  considered  in  this  study. 


2.4.  Analysis  of  laminar  flame  response 


In  the  following  section,  a  one-dimensional  analysis  is  pre¬ 
sented  to  investigate  the  response  of  the  flame  structure  to  dilu¬ 
tion.  This  analysis  provides  a  fundamental  understanding  about 
the  importance  of  considering  dilution  effects  on  the  combustion 
process. 


Fig.  6.  Steady  flamelet  solution  at  Z  =  Zc  for  different  dilution  levels  a  showing 
temperature  as  function  of  scalar  dissipation  rate.  The  operating  conditions 
correspond  to  the  MILD  combustion  experiment  of  [  1 6],  Fuel  is  methane  at  ambient 
temperature  and  air  is  preheated  to  600 1<.  Global  equivalence  ratio  is  equal  to 
0G  =  0.588  and  no  heat  losses  are  considered. 


2.4.1.  Dilution  without  heat-loss  effects 

Figure  6  shows  simulation  results  for  the  S-shaped  curve  with¬ 
out  consideration  of  heat  losses.  As  expected,  dilution  increases  the 
temperature  of  the  unburned  branch  and  the  temperature  of  the 
burnt  mixture  converges  to  a  unique  equilibrium  point.  The  critical 
conditions  for  the  quenching  and  ignition  scalar  dissipation  rates 
are  also  affected  by  dilution.  With  increasing  dilution  ratios,  the 
extinction  limit  increases  and  eventually  merges  with  the  ignition 
point.  The  transition  between  unburned  and  burnt  conditions  is 
then  smooth  and  monotonic,  which  is  in  accordance  with  the  MILD 
combustion  regime  characterization  of  Oberlack  &  Peters  [2],  Last, 
combustion  is  sustained  for  high  scalar  dissipation  rates  if  dilution 
is  sufficient:  this  explains  why  in  most  diluted  combustion  config¬ 
urations,  combustion  is  stabilized  despite  high  air  injection  veloc¬ 
ities.  The  influence  of  dilution  on  combustion  is  also  highlighted  by 
the  species  trajectories  in  compositional  space.  Figure  7  shows  the 
evolution  of  mass  fractions  of  carbon  dioxide  and  carbon  monoxide 
as  functions  of  Yc  for  different  dilution  rates  and  at  a  given  mixture 
fraction. 

From  these  results  it  can  be  seen  that  the  chemical  trajectories 
converge  to  an  attractor  space  that  is  defined  by  the  undiluted 
mixture  trajectory.  As  noted  previously,  the  equilibrium  composi¬ 
tion  is  identical  for  all  cases,  since  the  enthalpy  of  the  mixture  is 
not  affected  by  dilution.  However,  despite  the  convergence  to  the 


(a) 


same  manifold,  these  trajectories  are  strongly  dependent  on  the 
dilution  rate.  The  species  production  rates  are  also  strongly  depen¬ 
dent  on  the  dilution.  This  is  illustrated  in  Fig.  8b,  where  the  pro¬ 
gress  variable  and  the  CO  production  rate  are  plotted  as 
functions  of  Yc. 

2.4.2.  Dilution  with  heat-loss  effects 

In  practical  burner  configurations,  heat  losses  play  a  crucial  role 
in  the  combustion  process.  Figure  9  presents  various  sets  of  S- 
curves  obtained  for  different  values  of  a  and  diluent  temperatures 
TD|1.  Independent  of  the  level  of  heat  losses  applied  to  the  diluent, 
the  temperature  of  fresh  gases  increases  with  dilution.  The  equilib¬ 
rium  temperature,  obtained  for  low  scalar  dissipation  rates,  de¬ 
creases  when  TDl1  decreases  and  for  increasing  a.  For  each  case, 
the  unstable  branch  of  the  S-curve  progressively  vanishes  with 
dilution. 

The  scalar  dissipation  rate  value  under  quenching  conditions, 
Xqz,  changes  according  to  a  and  p.  Without  heat  losses  (case  (a)  in 
Fig.  9),  combustion  is  possible  at  very  high  strain  rates  if  dilution 
is  sufficiently  high:  the  increase  in  reactant  temperature  induced 
by  mixing  with  high-enthalpy  diluent  allows  chemical  reactions 
to  occur  even  under  high-strain-rate  conditions.  However,  for  dil¬ 
uent  temperatures  below  the  adiabatic  condition,  two  different 
scenarios  can  be  identified: 


(b) 


Fig.  7.  Profiles  of  species  mass  fractions  of  (a)  C02  and  (b)  CO  as  functions  of  Yc;  conditions  are  identical  to  those  used  in  Fig.  6. 


Temperature  [K]  Temperature  [K] 
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(a)  (b) 


Fig.  8.  Profiles  of  species  source  terms  of  (a)  Yc  and  (b)  CO  as  function  of  Yc;  conditions  are  identical  to  those  used  in  Fig.  6. 


(b) 


(c)  (d) 

Fig.  9.  Steady  flamelet  temperatures  for  Z  =  ZG  as  function  of  the  scalar  dissipation  rate  for  various  dilution  levels,  a.  Levels  of  increasing  dilution  are  indicated  by  the  arrow. 
Global  equivalence  ratio  is  equal  to  0.77.  (a)  Tm  =  T™0  =  2190  K;  (b)  Tm  =  T™02  =  1987  1<;  (c)  TDil  =  T™04  =  1766  K;  (d)  Tm  =  T™06  =  1530  K. 
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•  For  low-heat-loss  conditions,  the  trends  are  similar  to  the  adia¬ 
batic  dilution  case;  Xz  increases  with  dilution,  which  favors 
combustion  by  transferring  thermal  energy  to  the  system. 

•  If  heat  losses  exceed  a  certain  threshold,  Xz  decreases  with 
increasing  dilution  level.  As  result,  the  stable  branch  decreases 
accordingly. 

In  conclusion,  the  establishment  of  a  MILD  combustion  regime 
is  dependent  not  only  on  the  dilution  with  burned  reaction  prod¬ 
ucts  but  also  on  heat  losses.  While  a  pronounced  dilution  of  reac¬ 
tants  is  necessary  to  achieve  MILD  combustion,  heat  losses  are  of 
paramount  importance  in  affecting  the  flame  behavior  and  in 
extending  the  stability  range  in  confined  combustors  with  internal 
recirculation. 

3.  Turbulent  combustion  model 

3.3.  Presumed  probability  density  function  model 

A  statistical  description  is  used  to  represent  the  turbulence/ 
chemistry  interactions  at  the  sub-grid-scale  level.  A  presumed 
probability  density  function  (PDF)  approach  is  considered,  in  which 
the  filtered  scalar  quantities  are  written  as 

i }  =  Jjjj  G^,(Z,A,a,ll)P(Z,A,oc,fi)dZdAdadfS  (28) 

and  the  density-weighted  joint  PDF  P  is  expressed  as 

P(Z,  A,  a.,  ft)  =  =P(Z,  A,  a,  fl),  (29) 

with  P  being  the  nonweighted  joint  probability  density  function. 
The  progress  parameter  A  is  defined  to  be  statistically  independent 
of  the  mixture  fraction  Z  [31  ].  By  assuming  statistical  independence 
of  a  and  /3,  the  joint  PDF  can  then  be  expressed  in  terms  of  marginal 
distributions: 

P(Z,  A,  a, /?)  =  P(Z)P(A)P(a)P(/l).  (30) 


Transport  equations  for  the  first  two  moments  of  mixture  frac¬ 
tion,  progress  variable,  and  enthalpy  are  written  as 


^  +  V  •  (puZ)  =  V  .  (p(D  +  D  ~)  vz) , 
df>^  +V  ■{piiZZ2)  =  V  [p(D  +  D  ~)VZ^j 

(32) 

+  2pDT  \  VZ\2  -  sff, 

(33) 

+ V  •  (piiYc)  =  V  •  (p(D  +  Dt~  )  V?c)  +  -pdoYc, 

(34) 

+  V  •  (piiYd)  =  V  ■  (p(D  +  D~  )VY„)  +  pYod, 

(35) 

^  +  V  •  (puH)  =  V  •  (p(D  +  D-)  VH) . 

(36) 

The  mixture-fraction  subgrid  scalar  dissipation  rate  is  modeled 
using  a  linear  relaxation  assumption 


oV/2 

%=2  pCJh-p,  (37) 

in  which  A  is  the  LES  filter  size  and  Cx  is  a  model  constant  taken  as 
unity.  Radiation  effects  in  the  enthalpy  equation  are  neglected  and 
only  convective  heat  losses  to  the  combustor  wall  are  considered. 

To  solve  the  balance  equation  for  the  dilution  variable  (35),  the 
source  term  cod  needs  to  be  modeled.  This  term  accounts  for  inter¬ 
nal  exchange  of  reaction  products  from  the  compositional  stream 
to  the  dilution  stream.  Since  this  requires  the  consideration  of  mul¬ 
tispecies  interaction,  this  term  cannot  be  analytically  related. 
Therefore,  an  empirical  model  for  this  transfer  term  is  proposed. 
The  role  of  the  source  term  is  to  transform  burnt  gases  into  diluent. 
The  source  term  peaks  at  the  condition  corresponding  to  the  global 
equivalence  ratio.  Whether  or  not  the  incoming  streams  will  be¬ 
come  diluted  then  depends  on  the  convective  velocity  and  the  dif- 
fusivity.  For  this,  we  assume  that  burnt  products  are 
instantaneously  transformed  into  the  diluent  stream  in  the  burnt 
gas  regions.  Mathematically,  this  term  is  then  written  as 


In  the  following,  a  beta-PDF  is  used  to  model  the  mixture  frac¬ 
tion  distribution,  implemented  following  the  approach  of  Lien  et  al. 
[32  .  In  MILD  or  highly  diluted  combustion  regimes,  the  reaction 
zone  is  spatially  distributed,  which  is  the  reason  why  it  was  often 
compared  to  well-stirred  reactors  [33,34].  As  such,  the  distribu¬ 
tions  of  the  reaction  progress,  dilution,  and  heat-loss  variables 
are  represented  by  Dirac  functions.  After  A  is  replaced  by  the  pro¬ 
gress  variable  Yc,  the  following  thermochemical  state-space  repre¬ 
sentation  is  obtained  for  the  DFPV  model: 

xv(z,sz,yc,  a,/)).  (3i) 

Eq.  (31 )  is  used  to  provide  information  about  all  Favre-filtered  ther- 
mo-viscous-chemical  quantities  (for  example,  temperature,  species, 
source  terms,  density,  and  transport  properties).  It  is  noted  that  by 
introducing  Dirac  distributions  for  the  dilution  and  the  heat  loss 
parameters,  we  can  simplify  the  notation  and  now  write  a  =  a 
and  jl  =  [i. 

3.2.  Filtered  transport  equations 

In  addition  to  the  solution  of  the  conservation  equations  for 
mass  and  momentum,  five  additional  transport  equations  are  re¬ 
quired  to  close  the  system  of  equations  in  the  DFPV  model.  Here, 
scalar  subgrid  fluxes  are  modeled  by  introducing  turbulent  diffu- 
sivities  DT,h  which  are  evaluated  using  a  procedure  for  the  dy¬ 
namic  evaluation  of  the  turbulent  Schmidt  numbers  following 
the  method  of  Lilly  [35], 


H(C- 1).  (38) 

IV 

m 

This  expression  consists  of  four  terms  that  represent  the  follow¬ 
ing  physical  processes: 

•  Term  1:  This  term  is  the  inverse  of  the  time  needed  to  produce 
diluents,  which  is  here  evaluated  as  i=  At,  assuming  instanta¬ 
neous  transfer  of  products  into  the  dilution  stream. 

•  Term  II:  This  term  relaxes  YA  toward  its  equilibrium  value,  Y^'1 
is  dependent  on  Z  and  corresponds  to  the  value  of  the  product 
gas  composition  at  equilibrium.  Since  Yd  is  linearly  dependent 
on  Z,  it  can  be  shown  that  Y°l1  is  given  by: 

f0rS<^  (39, 

l''*ti-Sf',«U  for  z  >  zc. 

•  Term  III:  The  assumptions  made  earlier  state  that  the  diluent 
mixture  is  composed  of  burnt  gases  that  correspond  to  the 
product  mixture  at  the  global  equivalence  ratio.  The  Gaussian 
function  included  here,  centered  at  ZG,  localizes  the  production 
of  diluent  in  mixture  fraction  space  and  the  value  for  its  stan¬ 
dard  deviation  was  chosen  to  be  equal  to  a  =  0.05. 

•  Term  IV:  H(C  -  1)  is  the  Heaviside  function  on  C.  The  introduc¬ 
tion  of  this  function  makes  it  possible  to  form  diluents  only  in 
the  fully  burned  gas  mixture. 
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With  the  solution  of  Eqs.  (32)-(34),  the  DFPV-state-variables 
can  then  be  evaluated  as  follows: 


a  = 


Z  —  ocZg 
1  -  a 


Sz  = 
P  = 


z"2 

5 1  ~z)~ 

H  -  (Hi o  - 


Z"2 


(Z  -  Zc)(1  —  <X+  ccZc  -  Z) 
H°l0)(Z  -  ocZq)  -  (1  -  a)H°!0 


(40) 

(41) 

(42) 

(43) 


4.  Application  to  adiabatic  MILD  operating  regime 

To  examine  the  performances  of  the  DFPV  model  in  application 
to  an  adiabatic  MILD  combustion  regime,  large  eddy  simulations  of 
a  well-insulated  combustion  chamber  are  performed.  The  experi¬ 
mental  configuration  under  consideration  has  been  investigated 
by  Castela  et  al.  [25],  The  choice  of  using  LES  was  motivated  by  pre¬ 
vious  investigations  by  Graqa  et  al.  [36],  in  which  it  was  concluded 
that  RANS  computations,  associated  with  EDC  or  transported  PDF 
combustion  models,  provide  an  inadequate  description  of  the  sca¬ 
lar  mixing  processes  up  to  half  the  length  of  the  combustion 
chamber. 

In  this  investigation,  we  will  compare  the  following  three  com¬ 
bustion  models: 

•  The  laminar  flamelet  model  (LFM)  of  Peters  [18],  for  which 
only  the  stable  branch  of  the  S-curve  is  tabulated,  and 
extinction  due  to  quenching  is  not  captured. 

•  The  FPV  approach  [22],  which  considers  the  solution  of  the 
unstable  branch,  but  does  not  account  for  variations  in  the 
reactant  composition. 

•  The  DFPV  model  without  heat  losses  (p  =  0)  proposed  here. 


4.1.  Experimental  configuration 


diameter  of  90  mm.  A  schematic  of  this  configuration  is  shown  in 

Fig.  10. 

The  operating  conditions  investigated  here  are  summarized  in 
Table  1.  The  adiabatic  operating  condition  was  indirectly  con¬ 
firmed  by  comparing  computed  adiabatic  flame-temperature  re¬ 
sults  (at  the  global  equivalence  ratio)  with  reported  exhaust 
temperature  measurements,  indicating  that  the  measurements 
are  within  5%  of  the  theoretical  results.  We  thus  considered  that 
the  configuration  was  adiabatic. 


4.2.  Computational  setup 

The  finite-volume  code  YALES2  [38,39]  was  used  to  simulate 
the  combustor  configuration.  In  this  solver,  the  Navier-Stokes 
equations  are  solved  under  a  low-Mach  approximation  using  a  pro¬ 
jection  method  [40]  for  variable  density  flows.  For  the  presented 
computations,  the  equations  are  solved  in  conservative  form.  A 
fourth-order  finite-volume  scheme  is  used  to  discretize  the  spatial 
operators  and  a  fourth-order  Runge-Kutta-type  scheme  is  used  for 
time  advancement. 

To  investigate  the  sensitivity  to  the  LES  mesh  resolution,  simu¬ 
lations  on  three  different  computational  grids  were  performed,  and 
the  mesh  characteristics  are  summarized  in  Table  2. 

A  detailed  analysis  of  mean-flow  results  showed  that  species 
and  temperature  statistics  were  not  further  improved  when  the 
mesh  resolution  was  increased  from  the  medium  to  the  fine  grid 
(Fig.  If).  Therefore,  only  results  from  the  medium  mesh  are 
presented. 

The  lookup  table  was  discretized  with  201  x  15  x  201  x  11  grid 
points  in  the  direction  Z  x  Sz  x  C  x  a.  The  generation  of  the  filtered 
chemistry  table  was  achieved  in  approximately  2h  on  12 
processors. 

Inflow  conditions  for  the  oxidizer  stream  were  prescribed  from 
a  turbulent  velocity  profile  with  5%  turbulent  intensity.  An  analysis 
of  sensitivity  to  the  turbulent  intensity  was  carried  out,  and  negli¬ 
gible  impact  on  the  flow  field  was  observed.  As  for  the  progress 
variable,  the  boundary  conditions  values  associated  with  the  dilu¬ 
tion  variable  are  set  to  zero  for  every  inlet. 


The  configuration  is  a  reversed-flow  combustion  chamber,  in 
which  inlet  and  exhaust  ports  are  located  on  the  same  side.  This 
configuration  ensures  sufficiently  large  residence  times  to  com¬ 
plete  combustion  and  to  promote  intense  mixing  of  burned  gases 
with  the  unburned  reactant  stream  [37],  Similar  reverse-flow  con¬ 
figurations  have  previously  been  considered  [11,12]  in  order  to 
achieve  MILD  combustion  conditions. 

The  configuration  consists  of  a  cylindrical  combustion  chamber 
of  340  mm  length  and  100  mm  diameter.  The  injection  system  is 
composed  of  a  central  fuel  nozzle  with  a  diameter  of  4  mm,  supply¬ 
ing  natural  gas  in  the  volumetric  composition  ratio 
CH4/C2H6/C3H8/N2  =  85.1/7.6/1.9/5.4%.  Preheated  air  at  600  K 
is  supplied  through  an  annualar  injector  with  an  inner  diameter 
of  14  mm  and  an  outer  diameter  of  18.5  mm.  The  exhaust  port  is 
an  annulus,  having  an  inner  diameter  of  75  mm  and  an  outer 


4.3.  Instantaneous  flow  field 

A  comparison  of  instantaneous  temperature  field,  obtained 
from  all  three  models  is  shown  in  Fig.  12.  These  qualitative  results 
show  that  the  instantaneous  temperature  field,  predicted  by  LFM, 
shows  considerably  higher  temperatures  than  obtained  from  the 
other  two  models.  The  maximum  temperature  is  close  to  the  adia¬ 
batic  temperature  of  2150  K.  This  is  due  to  the  model  formulation, 
which  does  not  consider  flame  extinction  (only  the  upper  stable 
branch  of  the  S-curve  is  considered).  This  is  further  illustrated  by 
showing  scatter  plots  of  temperature  on  a  longitudinal  cut  in 
Fig.  13  (left).  Since  this  model  does  not  account  for  extinction 
events,  all  reactants  mix  and  burn,  and  eventually  relax  toward 
the  equilibrium  composition.  Combustion  is  therefore  considered 
to  be  fast. 


Fig.  10.  Schematic  view  of  the  combustion  chamber  geometry  of  Gra^a  et  al.  [36],  showing  (left)  a  longitudinal  cut  and  (right)  a  combustor  cross-sectional  view. 
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Table  1 

Global  equivalence  ratio  i pc,  inlet  velocity  of  air-  and  fuel  streams,  Uair  and  U fuei, 
respectively. 


<£g 

Uair  (m/s) 

Uftrel  (m/s) 

Reair 

Refuel 

Tg*  (S) 

0.416 

108 

17,7 

8751 

4514 

0.211 

Note:  Reynolds  number  based  on  air  injector  Reair  and  fuel  injector  RefllH,  and  global 
residence  time  Tgs  evaluated  as  the  ratio  of  the  volume  of  the  chamber  to  the  total 
reactant  volume  flow  rate. 


Table  2 

Mesh  information  for  the  simulation  of  the  configuration  of  Castela  et  al.  25). 


Mesh 

A0  (mm) 

DoF 

Elements 

Coarse 

0.4 

1.2  M 

6M 

Medium 

0.22 

5.5  M 

28.5  M 

Fine 

0.11 

41M 

228 M 

Note:  A0  denotes  the  minimum  filter  size  and  DoF  is  degrees  of  freedom  (in 
millions). 


The  FPV  approach  allows  for  the  presence  of  unburned  mixtures 
as  solutions  of  the  flamelet  equations  (Fig.  12  (middle)).  In  the 
present  case,  the  reactant  gas  temperature  at  the  inlet  is  too  low 
to  promote  autoignition.  Partial  premixing  of  reactants  occurs  after 
injection  into  the  combustion  chamber.  Reactants  are  then  diluted 
by  recirculating  hot  burned  gases,  so  that  the  flame  stabilization 
and  combustion  are  facilitated  by  diffusion  transfer  of  heat  and 


radical  species  to  the  partially  premixed  reactant  mixture.  The 
FPV  model  predicts  a  maximum  temperature  that  is  considerably 
lower  than  the  LFM  model.  This  is  further  illustrated  by  scatter 
plots  in  Fig.  13  (middle).  Because  of  the  partial  premixing  of  the 
reactants,  occurring  prior  to  combustion,  chemical  trajectories 
converge  to  the  equilibrium  condition  at  Z  =  Z^.,  thereby  bypass¬ 
ing  the  high-temperature  regions  around  the  stoichiometric  condi¬ 
tion.  It  is  noteworthy  to  point  out  that  the  FPV  model  predicts 
combustion  of  fuel-rich  pockets,  resulting  in  the  formation  of  hot 
spots  with  temperatures  around  1550  K;  these  spots  are  shown 
by  the  red  areas  in  Fig.  12  (middle). 

An  instantaneous  temperature  field  predicted  by  the  DFPV 
model  is  illustrated  in  Fig.  12  (right),  exhibiting  qualitative  similar¬ 
ities  to  the  temperature  field  that  is  predicted  by  the  FPV  model 
(Fig.  12  (middle)).  The  detachment  of  the  reaction  zone  from  the 
region  between  air  and  fuel  inlets  and  temperature  field  is  consid¬ 
erably  lower  than  that  obtained  from  the  other  two  simulations. 
Figure  13  (right)  shows  that  for  the  whole  mixture  fraction  range, 
unburned  gas  states  exist.  Moreover,  the  penetration  length  of  the 
reactants  is  lower  in  the  DFPV  model  than  in  the  FPV  model,  which 
highlights  the  effect  of  dilution  on  the  combustion  process. 

4.4.  Statistical  flow  field  results 

A  comparison  of  time-averaged  temperature  profiles  along  the 
combustor  centerline  is  shown  in  Fig.  14.  The  temperature  pre¬ 
dicted  by  the  LFM  model  is  overestimated  until  170  mm.  The  FPV 


Fig.  11.  Mean  temperature  and  dry  mole  fractions  of  C02  and  02  along  the  centerline  of  the  configuration  for  the  medium  and  fine  meshes. 


LFM  FPV  DFPV 


Fig.  12.  Temperature  fields  in  a  median  plane  for  the  three  models.  The  distribution  is  shown  close  to  the  injection  system  and  the  spatial  scale  is  given  in  meters.  Left:  LFM 
model;  Middle:  FPV  model;  Right:  DFPV  model. 
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Fig.  14.  Time-averaged  temperature  along  the  centerline  of  the  configuration. 
Symbols :  experimental  data  points;  Continuous  line:  DFPV  model;  Dashed  line:  FPV 
model;  Dashed-dot-dot  line:  LFM  model. 

model  predicts  a  lower  heat  release  in  the  nozzle-near  region,  and 
good  agreement  with  experimental  data  is  obtained  for 
x>150mm,  corresponding  to  44%  of  the  total  length  of  the 
combustion  chamber.  A  noticeable  improvement  is  obtained  with 
the  DFPV  model,  where  the  temperature  homogeneity  is  better 
reproduced.  It  is  still  unclear  why  discrepancies  arise  in  the  noz¬ 
zle-near  region,  but  it  can  be  speculated  that  the  mixing  process 
requires  more  precision,  because  of  measurement  uncertainties. 
Moreover,  it  is  noteworthy  that  the  discrepancies  appearing  in 
the  nozzle-near  region  are  significantly  lower  than  in  the  RANS 
computations  of  Grai^a  et  al.  [36],  where  the  temperature  at  the 
centerline  at  x  =  70  mm  is  overestimated  by  800  K. 

Similar  trends  are  observed  for  dry  mole  fractions  of  02  and  C02 
(shown  in  Fig.  15).  Considerable  discrepancies  can  be  observed  for 


Fig.  15.  Comparisons  of  temporally  averaged  centerline  profiles  for  mole  fractions 
of  02  and  dry  C02.  Symbols:  experimental  data  points;  Continuous  line:  DFPV  model; 
Dashed  line:  FPV  model;  Dashed-dot-dot  line:  LFM  model. 


the  predictions  with  the  LFM  model,  while  good  agreement  is  ob¬ 
tained  for  the  simulations  using  the  FPV  and  the  DFPV  models.  This 
provides  further  support  that  the  predicted  flame  characteristics 
near  the  injector  are  representative  for  this  burner  configuration. 

Comparisons  of  carbon  monoxide  mole  fractions  along  the  cen¬ 
terline  are  shown  in  Fig.  16.  Large  differences  are  observed  be¬ 
tween  the  three  model  predictions.  The  LFM  model  overpredicts 
the  reaction  progress  for  all  equivalence  ratios,  which  leads  to  a  ra¬ 
pid  increase  of  CO  in  the  nozzle-near  region.  This  is  followed  by  a 
rapid  consumption  of  CO  due  to  the  relaxation  towards  equilib¬ 
rium.  The  delayed  combustion  process  predicted  by  the  FPV  model 
yields  an  overprediction  of  CO.  In  contrast,  when  the  dilution  of 
reactants  is  considered,  the  model  qualitatively  captures  the  trend 
and  the  magnitude  of  the  CO  profiles. 

The  accurate  prediction  of  CO  is  particularly  challenging,  which 
is  due  to  its  sensitivity  to  the  reaction  chemistry  and  local 
flow-field  composition.  The  following  comparison  of  CO  results 
emphasizes  the  importance  of  considering  dilution  effects  in  the 
combustion  model.  To  illustrate  the  formation  and  consumption 
of  CO  inside  the  combustor,  we  evaluate  scatter  data  in  C0-C02 
composition  space  throughout  the  combustor.  These  results  are 
presented  in  Fig.  17,  confirming  that  dilution  has  a  considerable 
influence  on  the  CO  conversion.  When  dilution  effects  are 
neglected,  the  only  trajectory  available  is  defined  by  the  undiluted 
case,  and  the  production  of  CO  is  unsatisfactory.  The  results 
obtained  from  this  large  eddy  simulation  substantiate  the  laminar 
flame  investigations  and  show  that  the  flame  structure  and  minor- 
species  conversion  exhibit  pronounced  sensitivity  to  the  dilution  of 
the  reactant  gases.  Furthermore,  these  simulations  demonstrate 
the  capability  of  the  proposed  DFPV  model  in  predicting  internal 
dilution  systems  and  dilution-controlled  MILD  combustion. 


Fig.  16.  Comparisons  of  temporally  averaged  centerline  profile  for  dry  mole 
fractions  of  CO.  Symbols:  experimental  data  points;  Continuous  line:  DFPV  model; 
Dashed  line:  FPV  model;  Dashed-dot-dot  line:  LFM  model. 
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Fig.  17.  Scatterplots  of  CO  as  a  function  of  C02  (in  dry  mole  fractions)  for 
Z  =  ZG  ±0.002.  Black :  DFPV;  Red:  FPV;  Blue:  LFM. 


4.5.  Flame  structures  analysis 

The  Takeno  index  [41]  is  used  to  discriminate  combustion 
modes  in  the  combustion  chamber.  Yet  the  formulation  of  this 
index  poses  issues  when  detailed  chemistry  effects  are  taken  into 
account  and  in  the  absence  of  reactants.  To  overcome  these 
deficiencies,  another  flame  index  is  proposed. 

As  mentioned  in  several  studies  [42-44],  in  moderately  curved 
diffusion  flames,  the  gradients  of  2  (nz  =  VZ/|VZ|)  and  Yc 
( nYc  =  VYC/  |  VYC  |)  are  aligned  so  that  nz  ■  nYc  ±  1.  For  weakly 
stratified  flames  and  premixed  flames  that  are  subjected  to  low 
strain,  the  flame  front  is  perpendicular  to  iso  surfaces  of  mixture 
fraction,  which  implies  that  nz  ■  nYc  =  0.  In  partially  premixed 
combustion,  one  can  then  obtain  -1  <  nz  ■  hYc  <  1.  In  the  follow¬ 
ing,  we  define  the  flame  index  GIZ  as: 


VZ  ■  VYC 
“l  VZ  II  VYC| 


(44) 


To  ensure  that  we  analyze  only  the  reaction  zone,  this  index  is 
conditioned  on  the  source  term  of  the  reaction  progress  variable 
and  e  is  a  parameter,  which  is  here  chosen  to  be  equal  to  1  s  '. 
GIZ  =  ±1  defines  a  diffusion-controlled  combustion  mode.  It  is 
important  to  note  that  Eq.  (44)  does  not  account  for  subgrid  contri¬ 
butions.  To  assess  the  significance  of  subgrid  contributions,  we 
evaluated  the  flame  index  for  the  medium  and  the  fine  mesh, 
and  both  provide  quantitatively  similar  results.  The  flame  index 
is  computed  from  resolved  quantities.  It  is  expected  that  the  qual¬ 
itative  nature  of  the  predicted  flame  structure  is  weakly  sensitive 
to  the  errors  induced  by  the  turbulent  combustion  model,  even  if 
a  feedback  on  the  resolved  quantities  is  expected.  On  tetrahedra- 
based  meshes  with  high  variations  in  the  cell  size,  the  distribution 
of  GIZ  needs  to  be  weighted  by  cell  volumes  to  estimate  the  volu¬ 
metric  contribution  of  the  different  combustion  modes  that  are 
present  inside  the  burner.  The  discrete  volume-weighted  flame  in¬ 
dex  is  written  as 


EVa 


(z  —  ‘ 


cbc>e;GIz 


EE, 


(45) 


The  distribution  of  the  volume-weighted  flame  index  is  illus¬ 
trated  in  Fig.  18.  The  distribution  exhibits  a  bimodal  shape,  indicat¬ 
ing  that  combustion  mainly  occurs  in  the  diffusion-controlled 
regime.  A  quantitative  analysis  showed  that  more  than  55%  of 
the  reaction  zone  volume  is  characterized  by  Z-Yc  alignment  an¬ 
gles  in  the  range  of  [ — 7t/6,  7t/6]  or  [— 57t/6, 57t/6]  radians. 


GIZ  [-] 


Fig.  18.  Percentage  of  volume  occupied  by  cells  conditioned  by  the  flame  index  Glz. 


These  results  allow  us  to  conclude  that  the  combustion  in  this 
reverse-flow  combustor  is  primarily  dominated  by  nonpremixed 
combustion.  However,  a  nonnegligible  contribution  to  the  reaction 
occurs  in  gas  pockets  where  GIZ  differs  from  ±1,  indicating  the 
presence  of  premixed  and  partially  premixed  combustion  regimes. 
This  implies  that  the  DFPV  model  is  able  to  represent  these  flame 
topologies  at  strongly  diluted  MILD  combustion  regimes. 

5.  Application  to  a  nonadiabatic  MILD  operating  regime  under 
consideration  of  wall  heat-loss  effects 

In  the  previous  section,  the  applicability  of  the  DFPV  model  for 
the  prediction  of  dilution-controlled  combustion  regimes  without 
consideration  of  heat-loss  effects  was  examined.  However,  in  prac¬ 
tical  applications,  heat  losses  play  an  important  role  in  affecting 
the  combustion  stability  and  species  distribution.  The  perfor¬ 
mances  of  several  models  with  different  degree  of  complexity  are 
here  evaluated  in  application  to  a  nonadiabatic  configuration: 

•  The  FPV  approach  [22],  which  does  not  account  for  varia¬ 
tions  in  reactant  composition  and  heat-losses.  Note  that 
the  LFM  model  predictions  are  here  omitted  as  they  collapse 
on  the  FPV  approach  estimations. 

•  A  simplified  DFPV  approach,  denoted  “sDPFV,”  where  dilu¬ 
tion  and  heat  losses  are  correlated  in  time  and  space.  This 
approach  implies  a  unique  value  of  the  dilution  ratio  a  for 
a  given  enthalpy  defect.  This  model  is  derived  from  the 
DFPV  formulation,  and  modeling  details  are  presented  in 
Appendix  A. 

•  The  proposed  DFPV  model  allows  for  compositional  varia¬ 
tions  as  well  as  heat-loss  effects  without  enforced  correlations. 

5.1.  Experimental  configuration 

A  schematic  of  the  confined  combustor  configuration  is  pre¬ 
sented  in  Fig.  19.  This  burner  was  experimentally  studied  by  Veris- 
simo  et  al.  [16  .  The  cylindrical  combustor  has  a  diameter  of  50  mm 
and  a  length  of  50  mm.  The  reaction  products  exit  through  a  con¬ 
vergent  nozzle  of  length  150  mm  and  15  degrees.  Preheated  air 
at  673  K  is  injected  through  a  central  nozzle,  having  a  diameter 
of  10  mm.  Methane  fuel  under  ambient  conditions  is  supplied 
through  16  separated  nozzle  holes  (2  mm  in  diameter),  which  are 
coaxially  located  at  a  radius  of  15  mm.  The  fuel  inlet  velocity  is 
fixed  at  6.2  m/s,  and  the  experimental  investigations  consider  only 
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Fig.  19.  Schematic  view  of  the  combustion  chamber  geometry  of  Verissimo  et  al.  [16]  through  a  longitudinal  cut  (left),  and  a  cross-sectional  cut  through  the  burner  (right). 


Table  3 

Global  equivalence  ratio  0C,  air  inlet  velocities  Uair,  Reynolds  number  based  on  air 
injector  Re,,,,  and  global  residence  time  r]f\ 


Run 

0G 

d.ir  (m/s) 

Reair 

(S) 

2 

0.77 

113.2 

17,526 

0.29 

4 

0.58 

143 

22,140 

0.23 

variations  in  the  air  mass  flow  rate,  thereby  changing  the  overall 
equivalence  ratio. 

In  the  following,  two  operating  conditions  are  considered, 
namely  “run  2”  and  “run  4."  The  corresponding  operating  condi¬ 
tions  are  summarized  in  Table  3.  As  reported  in  the  experimental 
study,  high-momentum  air  injection  increases  the  recirculation 
and  the  mixing  of  unburned  reactants  with  reaction  products.  This 
results  in  a  homogeneous  temperature  field  and  ultralow  emis¬ 
sions  of  CO  and  NOx.  Verissimo  et  al.  [16]  studied  the  effect  of  ex¬ 
cess  air  on  the  combustion  regime  and  reported  a  gradual 
deviation  from  MILD  combustion  when  the  excess  air  coefficient 
was  increased.  Probe  measurements  were  conducted  to  determine 
species  concentrations  and  the  temperature  field  inside  the 
combustor.  From  the  measurements,  it  was  concluded  that  with 
increasing  excess  air  the  flame  structure  shifts  from  a 
homogeneous  reaction  zone  toward  a  recirculating  flame  that  sta¬ 
bilizes  in  proximity  to  the  fuel  injector. 

The  operating  point  denoted  as  “run  4”  features  the  character¬ 
istics  of  a  conventional  lean  combustion  process.  The  authors  of 
the  experimental  study  reported  a  localized  reaction  zone  and  a 
clearly  visible  flame  front.  For  the  second  case  (designated  "run 
2”)  the  measurement  indicates  that  the  combustion  occurs  in  the 
MILD  combustion  regime:  the  reaction  zone  is  distributed  over  a 
large  volume  and  no  flame  front  is  visible.  For  both  operating  con¬ 
ditions,  CO  and  NOx  emissions  are  below  10  ppm.  The  measure¬ 
ment  database  includes  mean  flow  results  for  temperature,  02, 
C02,  CO,  NOx,  and  unburned  hydrocarbons  (UHC)  along  axial  and 
radial  sections  in  the  burner.  The  reported  reproducibility  of  the 
experimental  data  is  5%  for  temperature  and  approximately  10% 
for  species  measurements. 

5.2.  Computational  setup 

Simulations  with  two  different  grid  resolutions  have  been  per¬ 
formed  (Table  4).  A  further  grid  refinement  study  was  performed 
but  did  not  show  noticeable  influence  on  the  scalar  mean  flow 
quantities.  Therefore,  only  results  obtained  with  the  fine  mesh  of 


Table  4 

Mesh  information  for  the  simulation  of  the  configuration  of  Verissimo  et  al.  [16], 


Mesh 

A0  (mm) 

DoF 

Elements 

Nair 

Coarse 

0.4 

2.58M 

13.7M 

-25 

Fine 

0.22 

8.45M 

48.9M 

-45 

Note:  A0  denotes  the  minimum  filter  size  and  DoF  is  degrees  of  freedom  (in  mil¬ 
lions);  Nllr  is  the  number  of  grid  points  along  a  diameter  of  the  air  inlet. 


48.9  million  control  volumes  are  presented  in  the  following.  The 
fuel  inlet  velocity  was  prescribed  by  a  laminar  profile  and 
turbulent  inlet  conditions  were  used  for  the  air  injection.  Turbulent 
fluctuations  with  different  turbulence  levels  (Tuinj  =  5%  for  run  4, 
and  Tuinj  =  5%  or  1%  for  run  2)  have  been  superposed  on  the  mean 
flow,  which  was  evaluated  from  a  periodic  pipe  flow  simulation. 
From  this  study  it  was  found  that  the  mean  velocity  profile  along 
the  centerline  is  sensible  to  the  turbulent  intensity,  but  no  exper¬ 
imental  information  about  the  velocity  field  was  reported  to  guide 
the  simulation  setup. 

Figure  20  shows  the  thermodynamic  equilibrium  temperature 
computed  for  different  values  of  the  equivalence  ratio.  The  range 
of  temperatures  accessed  in  the  experiments  is  added  in  the  figure 
for  the  </>G  =  0.58  and  </>G  =  0.77  cases.  The  maximal  temperature 
reached  in  the  experiment  is  150  K  and  350  I<  below  the  equilib¬ 
rium  states  for  run  4  and  run  2,  respectively.  This  comparison  im¬ 
plies  that  in  these  simulations,  the  impact  of  the  prescription  of  the 
boundary  conditions  is  of  importance,  since  heat  losses  are  signif¬ 
icant.  Wall  temperature  measurements  were  not  reported,  but  the 
extrapolation  of  the  measured  radial  temperature  profiles  from  the 
chamber  gives  an  estimate  of  the  wall  temperature  in  the  range  of 
1200  ±  200  K.  Based  on  this  estimate,  we  prescribed  a  constant 
wall  temperature  of  1200  K  for  all  simulations.  By  invoking  the 
low  temperature  of  the  combustion  process  and  the  relatively 
low  global  residence  times  found  in  the  configuration,  radiation  ef¬ 
fects  were  neglected  in  the  following  simulations. 

The  DFPV  table  was  discretized  with  151  x  15  x  121  x  11  x  11 
grid  points  in  the  direction  ZxSzxCxax/?.  The  generation  of 
the  DFPV  table,  including  the  computation  of  the  flamelets,  took 
approximately  5  h  on  24  processors. 


5.3.  Flow  field  topology 

The  high  momentum  of  the  air  stream  inlet  promotes  the  gen¬ 
eration  of  a  reverse  flow,  which  advects  hot  burnt  gases  toward  the 


4>H 


Fig.  20.  Thermochemical  equilibrium  temperature  as  a  function  of  the  mixture 
equivalence  ratio  $  under  Verissimo  et  al.  [16]  operating  conditions.  Arrows 
represent  the  range  of  measured  temperature  in  the  experimental  data,  and 
symbols  indicate  the  mean  temperature  at  the  exit  of  the  combustion  chamber. 
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fresh  reactants.  The  two  operating  conditions  impact  directly  the 
diluent  composition,  which  contains  large  amounts  of  oxygen  for 
run  4.  A  comparison  of  instantaneous  temperature  fields,  obtained 
under  the  two  operating  conditions,  is  shown  in  Fig.  21.  In  both 
cases,  low-velocity  regions,  localized  for  x  e  [0, 50]  mm  and 


Fig.  21.  Instantaneous  (left)  and  mean  (right)  temperature  fields  in  a  median  plane 
for  the  two  operating  conditions  using  the  DFPV  model.  Spatial  dimensions  are  in 
units  of  meters,  (a)  Run  4;  (b)  run  2. 


r  6  [20, 50]  mm,  exhibit  burned  gas  temperatures  lower  than  at 
the  chamber  exit.  This  is  attributed  to  the  heat  transfer  to  the  walls 
and  to  the  long  residence  times  found  in  these  regions.  For  run  4,  a 
high-mean-temperature  region  (>800  K)  can  be  located  at  the 
interface  between  the  oxidizer  stream  and  burned  gases  mixed 
with  methane.  For  run  2,  the  mean  temperature  distribution  is 
more  homogeneous  and  the  temperature  is  higher  than  in  run  4 
due  to  its  global  equivalence  ratio,  which  is  closer  to  stoichiometry. 
The  instantaneous  distributions  of  run  2  present  localized  and 
intermittent  high-temperature  spots  (>1900  I<).  The  mean  behav¬ 
ior  of  the  two  cases  is  similar  to  that  observed  experimentally. 

5.4.  Statistical  flow-field  results  for  run  4:  =  0.58 

Radial  profiles  of  temperature  and  dry  species  mole  fractions 
are  given  in  Fig.  22  for  run  4.  Experimental  results  are  shown  by 
symbols.  The  predicted  temperature,  obtained  from  the  two  DFPV 
approaches,  is  in  good  agreement  with  measurements.  The  FPV 
model  overestimates  the  temperature  by  400  K  throughout  the 
combustion  chamber.  This  shows  the  importance  of  heat-loss 
and  dilution  effects.  Temperature  predictions  that  are  obtained 
from  both  DFPV  models  exhibit  moderate  discrepancies  (<100  K) 
for  radial  position  r  >  5  mm.  By  computing  unsteady  flamelets 
subject  to  radiation,  comparison  between  the  characteristic  radia¬ 
tion  time  [31]  and  the  global  residence  time  of  the  configuration 
showed  that  temperature  overpredictions  can  be  attributed  to 
the  neglect  of  radiation  effects. 

Predictions  of  major  species  mole  fractions  are  in  good  agree¬ 
ment  for  all  models.  Contrary  to  the  temperature  distributions,  ma¬ 
jor  species  are  less  sensitive  to  heat  losses,  which  explains  the 
similar  results.  However,  discrepancies  can  be  observed  along  the 
centerline  of  the  configuration,  for  x  e  [45, 113]  mm  and  r  e  [0, 5] 
mm.  The  mole  fraction  of  02  is  underestimated,  while  C02  is 
overestimated. 

Similarly  to  the  results  presented  in  Section  4,  the  predictions  of 
carbon  monoxide  feature  very  different  evolutions  depending  on 
the  combustion  model.  In  the  experimental  data,  the  CO  spreads 
widely  in  radial  direction  at  the  first  few  probe  stations  (from  45 


Fig.  22.  Comparison  of  statistical  results  at  different  axial  locations  in  the  burner  for  the  operating  conditions  of  run  4,  showing  (from  top  to  bottom)  temperature,  and  dry 
mole  fractions  of  02,  C02  and  CO.  Symbols:  experimental  data  points;  continuous  line :  DFPV  model;  dashed-dot-dot  line :  sDFPV  model;  dashed  line:  FPV  model. 
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Fig.  23.  Comparison  of  statistical  results  at  different  axial  locations  in  the  burner  for  the  operating  conditions  of  run  2,  showing  (from  top  to  bottom)  temperature,  and  dry 
mole  fractions  of  02,  C02,  and  CO.  Symbols:  experimental  data  points;  continuous  line:  DFPV  model;  dashed-dot-dot  line:  sDFPV  model;  dashed  line:  FPV  model. 


to  147  mm).  The  DFPV  model  is  the  only  model  able  to  reproduce 
the  correct  order  of  magnitude  for  CO,  even  if  the  predictions  are 
not  satisfactory  from  79  to  147  mm.  Yet  Fig.  22  provides  evidence 
of  the  importance  of  the  product  gas  dissociation  and  heat  losses. 
More  specifically  for  pollutant  emissions,  the  observed  differences 
between  the  DFPV  models  can  be  attributed  to  the  prescription  of 
the  f)  parameter  in  the  sDFPV  model,  which  implies  that  several 
chemical  trajectories  are  not  considered  (see  A.2  for  further 
details). 

5.5.  Statistical  flow-field  results  for  run  2:  <^c  =  0.77 

Radial  profiles  of  temperature  and  dry  species  mole  fractions  for 
run  2,  corresponding  to  a  MILD  combustion  regime,  are  compared 
with  experimental  data  in  Fig.  23.  Similarly  to  run  4,  the  FPV  model 
without  consideration  of  dilution  and  heat  losses  overpredicts 
temperature  and  CO  profiles.  Temperature  predictions  from  the 
DFPV  and  the  sDFPV  models  are  very  similar,  and  in  good  agree¬ 
ment  with  experimental  data.  It  is  encouraging  to  note  that  both 
developed  models  converge  toward  the  same  equilibrium  (probe 
stations  located  at  xe[181  mm;  310  mm]).  However,  a  slight 
underestimation  of  temperature  is  noticeable  along  the  axis  of 
symmetry  for  regions  located  between  45  and  79  mm.  Away  from 
the  core  region  (r  >  5  mm),  the  DFPV  models  predict  higher  tem¬ 
peratures  than  seen  in  the  experiments.  Computations  of  unsteady 
laminar  flamelets  subject  to  radiation  for  these  operating  condi¬ 
tions  suggest  that  this  might  be  due  to  the  neglect  of  radiation. 
However,  these  discrepancies  might  also  be  sensitive  to  the  pre¬ 
scribed  wall  temperature,  for  which  no  experimental  data  are 
available.  Yet  these  models  demonstrate  their  capacity  to  repro¬ 
duce  heat  loss  effects. 

Predictions  for  major  species  are  in  better  agreement.  Discrep¬ 
ancies  arise  in  the  region  x  e  [11.113]  mm  and  re  [5,15]  mm,  for 
which  02  is  under-predicted  and  C02  is  overestimated.  A  similar 
discrepancy  between  temperature  and  major  species  was  also  ob¬ 
served  for  run  4.  These  discrepancies  may  be  attributed  to  the 
modeling  hypothesis  on  the  diluent  composition,  or  to  experimen¬ 
tal  uncertainties.  Predictions  of  CO  are  consistent  with 


temperature  and  major  species  estimate.  The  missprediction  of 
the  CO  mole  fraction  at  the  first  measurement  station  can  be 
explained  by  the  sensitivity  of  CO  to  the  temperature  field. 

Differences  between  results  from  the  DFPV  and  the  sDFPV  mod¬ 
els  are  very  low  for  these  operating  conditions.  For  run  2,  the  DFPV 
model  relaxes  toward  the  sDFPV  approach,  indicating  that  for  this 
specific  condition  the  combustion  process  is  most  likely  controlled 
by  heat  losses  rather  than  dilution  effects.  Nevertheless,  both  DFPV 
approaches  are  able  to  adequately  reproduce  the  first  moments  of 
species  and  temperature  field  for  both  conventional  and  MILD 
combustion  regimes. 

6.  Conclusions 

Modeling  challenges  of  MILD  combustion  arise  from  the  pro¬ 
nounced  dilution  of  reactants,  which  induces  compositional  varia¬ 
tions.  The  dilution  influences  the  reaction  chemistry  and  models 
should  therefore  be  able  to  account  for  these  effects.  Moreover, 
most  practical  combustor  configurations  operate  under  nonadia- 
batic  conditions,  which  requires  the  consideration  of  wall  heat-loss 
effects  to  adequately  describe  the  flame  stabilization,  fuel  conver¬ 
sion,  and  emissions. 

In  this  contribution,  a  turbulent  combustion  model  was  pro¬ 
posed  that  extends  a  tabulated  chemistry  model  to  account  for 
product  gas  dilution  and  heat-loss  effects.  Simplifying  assumptions 
were  formulated  to  decrease  the  model  complexity.  Two  tabula¬ 
tion  approaches  that  extend  the  FPV  model  were  introduced:  the 
DFPV  model,  where  dilution  and  heat  losses  may  evolve  on  differ¬ 
ent  spatiotemporal  scales,  and  the  sDFPV  model,  where  both  phe¬ 
nomena  are  constrained.  To  investigate  effects  of  dilution  and  heat 
losses  on  the  flame  structure,  a  laminar  flame-response  analysis 
was  performed.  Results  from  this  analysis  emphasize  the  signifi¬ 
cance  of  dilution  in  the  nonequilibrium  chemistry  and  flame 
stabilization. 

A  turbulent  closure  model  was  proposed,  using  a  presumed  PDF 
method  to  account  for  turbulence  chemistry  coupling.  An  exchange 
model  for  transferring  reaction  products  into  the  dilution  stream 
was  formulated.  This  model  was  then  applied  in  large  eddy 
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simulations,  and  two  different  combustor  configurations  were  con¬ 
sidered.  In  the  first  case,  results  obtained  with  the  DFPV  model 
were  in  good  agreement  with  experimental  data.  It  was  shown  that 
the  consideration  of  dilution  effects  leads  to  considerable  improve¬ 
ments  in  the  CO  prediction,  demonstrating  that  the  thermochemi¬ 
cal  trajectory  is  influenced  by  dilution.  In  the  second  application,  a 
combustor  configuration  was  considered  that  is  operated  in  the 
conventional  and  the  MILD  operating  regime.  Overall,  the  simula¬ 
tion  results  obtained  with  the  proposed  model  are  in  reasonable 
agreement,  demonstrating  the  ability  of  the  model  to  capture  ef¬ 
fects  of  species  dilution  and  heat  losses  on  the  flow-field  structure 
and  species  composition.  Future  improvements  in  the  model  re¬ 
quire  the  consideration  of  radiation  effects  and  the  generalization 
of  the  dilution  exchange,  which  is  currently  represented  by  an 
empirical  closure  model. 

The  proposed  tabulated  chemistry  method  assumes  that  the  in¬ 
ner  structure  of  the  turbulent  flame  front  is  close  to  a  laminar 
flame.  This  assumption  remains  to  be  further  investigated  for 
low-temperature  flames.  Indeed,  when  the  chemical  activities  slow 
down,  turbulent  eddies  penetrate  the  flame  front  and  may  affect 
the  chemical  trajectories.  Modeling  this  phenomenon  is  a  challeng¬ 
ing  problem. 
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Appendix  A.  Simplified  DFPV  approach  with  correlated  dilution 
and  heat-loss  evolution 


A.  3.  Model  formulation 


A  simplified  DFPV  tabulation  approach,  similar  to  the  one  devel¬ 
oped  by  Abtahizadeh  et  al.  [28],  can  be  obtained  by  assuming  that 
the  temperature  of  the  diluent  is  fixed  (/?  is  therefore  assumed  con¬ 
stant  in  the  entire  combustor).  A  schematic  of  this  approach  is 
illustrated  in  Fig.  A.l.  This  assumption  constrains  heat  losses  to 
the  dilution  rate  and  implies  that  both  phenomena  evolve  similarly 
in  space  and  time.  When  the  diluent  temperature  is  fixed,  the  reac¬ 
tion  chemistry  is  then  tabulated  as  a  function  of  three  parameters: 

^=^(Z,YC,  a).  (A.l) 

In  the  following,  this  model  is  referred  to  as  the  “sDFPV”  model. 
The  validity  of  this  simplifying  assumption  is  assessed  in  an 
application  to  a  practical  configurations  in  Section  5.  In  this  model, 
heat  losses  and  dilution  effects  are  linked,  so  that  H0'1  is  constant 
and  a.  reduces  to 


a  = 


H-TYLjF  ljOx  \  rrOx 

_ **<y.=o)  ^a=0 

HDil  o-  /  rjF  rrOx  \  r_rC>x 

^G(,na=0  Ha=0)  na=0 


(A.2) 


In  the  context  of  LES,  the  filtered  relation  for  a  is  directly 
expressed  from  Eq.  A.2, 


a 


7?  /  ljF  ljOx  \  rjOx 

_ **qc=q)  Hqi=0 

HDil  cr  ( LfF  rjOx  \  rrOx  5 

^G{na=o  n  a=0J  na=0 


(A.3) 


and  this  model  requires  only  information  about  the  enthalpy  H  and 
the  mixture  fraction  Z. 


a 


Fig.  A.l.  Schematic  representation  of  the  different  constraints  imposed  on  temperature  and  species  of  the  boundary  conditions  used  to  generate  counterflow  diffusion 
flames;  left :  DFPV  model;  right:  sDFPV  model. 


(a)  (b) 


Fig.  A.2.  Flamelet  solutions  for  Z  =  ZG  (<f>  =  0.58)  and  a  diluted  reactant  streams  temperature  of  1775  K  for  different  values  of  a  and  /?.  For  the  sake  of  clarity,  only  three 
solutions  have  been  illustrated,  (a)  Reaction  rate  of  progress  variable;  (b)  CO  mass  fraction. 
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A.2.  A  priori  comparison  between  DFPV  and  sDFPV  models 

It  is  noteworthy  that  the  sDFPV  approach  depends  on  the 
chosen  value  of  p.  A  preliminary  analysis  is  done  to  assess  the  sen¬ 
sitivity  of  the  flame  responses  when  different  values  of  the  heat 
loss  parameter  are  applied.  To  analyze  the  importance  of  neglected 
trajectories  in  the  sDFPV  model,  the  study  of  the  chemical  trajecto¬ 
ries  giving  the  same  diluted  reactant  streams  temperature  is  of 
interest.  By  prescribing  a  temperature,  at  a  given  equivalence  ratio, 
an  iterative  sorting  procedure  can  be  employed  to  recover  the  cor¬ 
responding  (a,  p)  parameter  pair. 

In  the  following,  the  mixture  fraction  is  fixed  to  Z  =  ZG  for  the 
run  4  investigated  in  Section  5.1,  and  the  diluted  reactant  streams 
temperature  is  fixed  to  1775  K,  which  corresponds  to  the  maximal 
mean  temperature  obtained  in  this  case.  Figure  A.2(a)  illustrates 
three  solutions  of  the  sorting  procedure  and  presents  the  reaction 
rate  6jYc  as  a  function  of  Yc.  The  DFPV  model  accounts  for  all  solu¬ 
tions,  whereas  the  sDPFV  model  allows  only  for  a  single  trajectory. 
The  strongest  reaction  zone  is  shifted  toward  high  values  of  Yc  for 
increasing  values  of  the  diluent  temperature.  As  these  trajectories 
are  strongly  dependent  on  the  parameterization  (a,  p),  the  reaction 
zone  can  consequently  be  poorly  predicted  by  the  sDFPV  model  if 
nonequilibrium  chemistry  effects  are  important. 

The  influence  of  the  diluent  temperature  on  species  evolutions 
is  highlighted  in  Fig.  A.2(b).  Even  if  the  equilibrium  mass  fractions 
are  very  similar,  very  different  trajectories  are  found.  Therefore,  the 
arbitrary  prescription  of  the  value  of  p  required  in  the  sDFPV  model 
can  have  a  dramatic  influence  on  the  combustion  process  and  pol¬ 
lutant  formation. 
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